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Maximum Entropy Method (MEM) is applied to glueball correlators at finite temperature constructed by using 
5,500 - 9,900 gauge field configurations generated in SU(3) quenched lattice QCD with the lattice parameter 
Aat = 2Nc/g = 6.25 and the renormalized anisotropy as/at — 4. The results support the thermal width 
broadening of the 0^"'' glueball near the critical temperature. 



It is believed that, in the neighborhood of the 
critical temperature T^, the QCD vacuum be- 
gins to change its structure such as the reduction 
of the string tension, the partial restoration of 
spontaneous chiral-symmetry breaking and so on. 
Since the hadrons are composites of quarks and 
gluons, whose interactions depend on the prop- 
erties of the QCD vacuum, some of the hadrons 
are expected to change their structure drastically 
near Tc- For instance, effective model studies sug- 
gest the pole-mass reductions of charmonium and 
light qq mesons as consequences of the reduction 
of the string tension and the partial restoration of 
spontaneous breaking of chiral symmetry, respec- 
tively |ll2l3j . In fact, these changes are consid- 
ered as important prccritical phenomena of QCD 
phase transition in RHIC QGP experiments, and 
corresponding lattice QCD calculations have been 
performed at quenched level |4I5) . In this paper, 
we study the thermal 0++ glueball, whose mass 
reduction seems natural in terms of the large dif- 
ference between the glueball mass tog — 1.5GeV 
and Tc ~ 280MeV in quenched QCD p7j . 

To study the pole mass of a hadron in lattice 
QCD, one has first to construct a temporal cor- 
relator as G{t) = (0(t)0(O)), and then resorts to 
its spectral representation as 



G(r) 



dLuK{T,uj)A{uj), 



where K{t,lu) 



cosh(cj(;3/2-T)) 
sinh(/3Lj/2) ' 



(1) 



/3 = 1/T, and 



A{uj) is the spectral function with its spatial mo- 
mentum projected to zero, i.e., A(w) = A{uj,p — 
0). Each peak position of A(lo) provides us with 
a pole mass of a hadron at T > 0, which can be 
observed in high-energy experiments. 

The most popular way to extract the spectral 
function A(lS) from the constructed correlator is 
to adopt an ansatz for A(lo), and then to per- 
form the fitting analysis. In the case of glue- 
balls, we first investigated their pole masses us- 
ing the simple narrow-peak ansatz with A{tj}) = 
C {S{lu — m) ~ S{uj + m)}, where C and to are 
the strength and the temperature-dependent pole 
mass, respectively 

We then proceeded to an advanced analysis 
adopting the Breit-Wigner ansatz, a more gen- 
eral and sophisticated ansatz, for the shape of 
the spectral function as [H] 

A{lj) = C {Sriio - LJo) - SriiJ + uJo)} (2) 



expressed with the peak function of 
dr{uj) = 



TT Cj2 + r2 ■ 



(3) 



Here, C, loq and T are the strength, the peak 
center corresponding to the pole mass, and the 
thermal width of the glueball peak, respectively. 
(The thermal width should not be confused with 
the decay width. It is important to keep in mind 
that, at r > 0, the thermal width is generated 
even for stable particles through the interaction 
with the thermally-excited particles.) 
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The Breit-Wigner analysis indicates that the 
main thermal effect for the 0''"+ glueball appears 
as the thermal width broadening [3] of about 
300 MeV with a reduction in the peak center of 
about 100 MeV. This Breit-Wigner ansatz for the 
lowest-lying peak in the spectral function seems 
natural and solid at relatively low temperature. 
However, near and above the critical tempera- 
ture, more complicated structure may appear in 
the spectral function. 

In this point, Maximum Entropy Method 
(MEM) is quite attractive, because it provides 
us with a numerical procedure to reconstruct the 
spectral function A{lu) directly from lattice QCD 
Monte Carlo calculations |10lll| . In this paper, 
we present our preliminary results on the MEM- 
reconstructed 0++ glueball spectrum at T > 0. 

We first consider the smearing method. Use 
of the extended operator is one of the most use- 
ful techniques to enhance the low-energy spec- 
tra. However, it has a disadvantage that it may 
make the physical interpretation of resulting peak 
less trivial. It can create an unphysical bump 
structure even in the spectral functions of non- 
interacting particles For instance, we con- 
sider the extreme case for glueballs. The naive 
continuum limit of the O"*""*" glueball operator is 
given as (j) oc C^jC^j. After the smearing, it is 
spatially extended according to the Gaussian dis- 
tribution as 



oc 



(27r)3/2p^ 



exp 



2p^ 



GUy)G1Az), 



where p controls the size of the spatial exten- 
sion P] . We adopt here the Coulomb gauge to ob- 
tain the continuum expression. We calculate the 
corresponding two-point functions perturbatively 
up to 0(a2), and extract their spectral functions. 
While the non-smeared spectral function is given 
as A{ui) (x w^, the smeared spectral function ac- 
quires an additional Gaussian factor as 



A{uj) oc uj^ exp {-{ujp)^/4} 



(5) 



which possesses an unphysical bump a,t lu — 
2y/2/p. The smearing method suppresses the 
overlap to the higher spectral components, and, 
as a consequence, this factor appears. 



The actual low-lying glueball is not a perturba- 
tive system but a definite bound state/resonance 
in quenched QCD below Tc- Note that, in order 
to study the mass and the width of a difinite res- 
onance, the problem of unphysical bump is not 
serious, because the pole position in the complex 
LU plane is unaffected by the smearing In this 
paper, we use the smearing method to enhance 
the lowest-lying 0++ glueball peak, which is neg- 
ligibly small in the non-smeared spectral function. 

We use the SU(3) anisotropic lattice plaquette 
action as 



Sc 



J2 ReTr{l-P,,(s)} 



Aat^ 

Nc 7G 

+ §r^G E ReTr{l-F,4(s)}, 



(6) 



s,t< j<3 



s.i,j<3 



where P^v{s) € SU(3) denotes the plaquette op- 
erator in the /i-i'-plane. The lattice parameter 
and the bare anisotropic parameter are fixed as 
Aat = 2iVc/5^ = 6.25 and 70 = 3.2552, re- 
spectively, so as to reproduce the renormalized 
anisotropy as Qs/at — 4. The scale unit is intro- 
duced from the on-axis data of the static inter- 
quark potential with the string tension ^/a — 440 
MeV. The resulting lattice spacings are given as 
= 9.365(66) GeV and aj^ = 2.341(16) GeV. 
The critical temperature is estimated as Tc ~ 280 
MeV from the susceptibility of the Polyakov loop. 
We generate 5,500-9,900 gauge configurations to 
(^4^ construct the glueball correlators, where statisti- 
cal data are divided into bins of the size 100 to 
reduce possible auto-correlations near the critical 
temperature. We adopt an appropriate smearing 
to enhance the low-energy spectrum. 

To reconstruct the spectral function A{uj), we 
apply MEM to the appropriately smeared glue- 
ball correlator normalized according to G{t = 
0) = 1. As the practical numerical procedure, 
we in principle follow Ref. ^I]. We adopt the 
following Shannon- Jaynes entropy as 



S: 



A{lu) - m(w) - A{lu) \oi 



\m{LL!) 



A7) 



where m{uj) is a real and positive function re- 
ferred to as the default model function. m(u;) 
is required to mimic the asymptotic behavior of 
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Figure 1. Reconstructed spectral functions of the 
lowest 0+"'' glueball from the temporal correlators 
atT= 130,253,275 MeV. 

A{uj) as w — > GO. As m{uj), we adopt Eq.(|SJ), i.e., 
the perturbative spectral function up to 0{a'^) as 

to(w) = iVw'*exp{-(wp)V4} , (8) 

where the normalization factor N is determined 
so as to mimic G(t = 0) = 1, i.e., 

POO 

1 = / dujK{T = 0,w)m(w). (9) 
Jo 

In Fig. n we show the reconstructed spec- 
tral functions of the lowest 0++ glueball from 
the temporal correlations at T = 130, 253, 275 
MeV. Since the error bar estimated by following 
Ref. jlll appears to be unreasonably small, we do 
not put it in Fig. ^ to avoid unnecessary confu- 
sion. For a reasonable estimate, it seems neces- 
sary to use the jackknife error estimate ^3]. In 
Fig.^ we see the tendency that the peak becomes 
broader with increasing temperature. 

To summarize, we have applied Maxi- 
mum Entropy Method (MEM) to appropriately 
smeared glueball correlators constructed with the 
anisotropic SU(3) lattice QCD at finite temper- 
ature below Tc, and have presented our prelimi- 
nary results. We have seen the tendency that the 
peak becomes broader with increasing tempera- 
ture, which supports the thermal width broaden- 



ing PI of the glueball near the critical tempera- 
ture. It is interesting to apply MEM to the glue- 
ball correlator above Tc, where it becomes less 
trivial to find a proper ansatz for the spectral 
function. Note that there may survive some non- 
perturbative effects even above the critical tem- 
perature. In fact, it was suggested, in a-n sector, 
that a strong correlation may still survive above 
the critical temperature ^3]. To investigate the 
correlation above Tc with MEM, it would be nec- 
essary to analyze the non-smeared correlators. 
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